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We present analysis of the cross-over behavior of disordered interacting two-dimensional electron systems 
in the parallel magnetic field. Using the so-called cross-over one-loop renormalization group equations for 
the resistance and electron-electron interaction amplitudes we qualitatively explain experimentally observed 
transformation of the temperature dependence of the resistance from a reentrant (nonmonotonic) behavior in 
relatively weak fields into an insulating-type behavior in stronger fields. 
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Disordered two-dimensional electron systems (2DES) 
have been in the focus of experimental and theoreti- 
cal research for several decades [1]. The interest to 
'2DES has been renewed because of the discovery of 
,the metal-insulator transition (MIT) in a high mobility 
'silicon metal-oxide-semiconductor field-effect transistor 
(Si-MOSFET) [2]. During last decade an interesting 
behavior of resistance and spin susceptibility has been 
found experimentally not only in Si-MOSFET but in 
pther 2D electron systems [3]. Recently, a major step 
'toward the theoretical proof for the MIT existence in 
2DES has been made in Ref. [4] . 

' If an electron density is higher than the critical one 
.(metallic phase) at low temperatures T ^ [r stands 
'for the transport mean-free path time] the increase of 
■the resistance with decreasing temperature is replaced 
|by the drop as T becomes lower than some Tmax (see 
'Fig. [1]) [3]. This nonmonotonic behavior of the re- 
.sistance has been predicted from the renormalization 
'group (RG) analysis of the interplay between disorder 
and electron-electron interaction in 2DES [5,6]. As a 
weak magnetic field _B|| is applied parallel to 2DES the 
'decrease of the resistance is stopped at some temper- 
mature and the resistance increases again [7]. Further 
increase of leads to the monotonic growth of the re- 
sistance as temperature is lowered. 

In the Letter we present the theoretical explanation 
for this striking behavior of the resistance in parallel 
magnetic field. We demonstrate that it can be ex- 
plained with the help of the RG analysis of disorder and 
electron-electron interaction in 2DES in the presence of 
the Zeeman splitting. 

The presence of the parallel magnetic field introduces 
a new energy scale gL ^-sB^ into the problem. Here gL 
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and /i B stand for the Lande factor and the Bohr magne- 
ton respectively. The Zeeman splitting g^^sBw sets the 
cut-off for a pole in the diffusive modes (diffusons) with 
opposite electron spin projections [5]. In the tempera- 
ture range ql^Ji-bBw <C T <C t^^ this cut-off is irrelevant 
and 2DEG behaves as if no parallel magnetic field is ap- 
plied. However, at lower temperatures T <^ g^^sB^^ the 
diffusive modes with opposite electron spin projections 
do not contribute and 2DES behaves as in the presence 
of the strong parallel magnetic field B^ ~ {gLHBT)~^ . 
It is well known [5] that in this case the resistance mono- 
tonically grows as T is lowered. Therefore, if the paral- 
lel magnetic field is weak qliibB^^ <^ Tmax! then the T- 
dependence of the resistance manifests the reentrant be- 
havior with maximum and minimum. At higher paral- 
lel magnetic fields ql^bB^ ^ Tmax the insulating-type 
behavior of the resistance restores such that it grows 
monotonically with decreasing T (see Fig. [1]). 

Below we shall illustrate this simple physical explana- 
tion of non-monotonic behavior of the resistance by the 
p 




Fig. [T] The sketch of the temperature dependence of 
the resistance in the presence of _B|| . The curve a) 
corresponds to the case _B|| = 0, the curve b) to the 
case gL/isBii <^ Tmax and the curve c) to the case 
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so-called cross-over RG equations (cf. Eqs. ([T7 1) - (PT|) ) 
for the resistance and the electron-electron interaction 
amplitudes that smoothly interpolate between the well- 
known cases of — and strong field. 

We consider two-dimensional interacting electrons 
with riy valleys in the presence of the quenched dis- 
order and the parallel magnetic field at low temper- 
atures T <C T^^. We assume that a magnetic field 
B± > T/ (De) where e and D stand for the electron 
charge and diffusion coefhcient respectively is applied 
perpendicular to 2DES in order to suppress the Cooper 
channel. We suppose that both the inverse intervalley 
scattering time and the valley splitting are much less 
than the temperature. 

Following Finkelstein [5], the effective quantum the- 
ory of disordered interacting two-dimensional electrons 
is given in terms of the generalized non-linear cr-model 
involving unitary matrix field variables Q^n'^'^^''^ 
which obey the constraint (r) = 1 . Here the integers 

= 1, 2, . . . , Nr denote the replica indices, m, n corre- 
pond to the discrete set of the Matsubara frequencies 
Un = ■KT{2n + 1). The integers (i = ±1, ±2, • • ■ , ±ny 
are the combined spin and valley indices. The effective 
action is given as 



S 



Here Ca is the free electron part 



tr(VQ)2 



(2) 



where a denotes the mean-field conductance in units 
/h with h being the Plank constant and symbol tr is 
the trace over Matsubara, replica, spin and valley in- 
dices. The Cf involves the electron-electron interaction 
amplitudes which describe the scattering on small (F) 
and large (r2) angles and the quantity z (originally in- 
troduced by Finkelstein [5]) which is responsible for the 
specific heat renormalization [9], 

Cf = 47rTztrr/(A- Q) -7rrF^tr/^Qtr/^„Q 

ail 

+7rTF2^(tr/^Q) ® (tr/^„Q) + ii^Tzirj^K. (3) 

an 

Here tr A(8)trS ^""'''"'"-S,™™''"''' and the matrices A, 
77 and II are given as 



(4) 



In the presence of the parallel magnetic field _B|| the 
Zeeman splitting should be taken into account [5] 



Cb, = -^Z29L^iBB\\ tr r,g + ^hd^^^sl (5) 

Here Z2 = z + T2 and the Pauli matrix is defined as 

{r,)Zi:'--i'^' =sign{Ci)dnmS''^Si^i\ (6) 

We mention that the last term in Cb^^ corresponds to 
the Fermi-liquid spin susceptibility. Finally, the term 



Ch 



4n,„ 



trAQ 
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is not a part of the theory but we shall use it later on 
as a convenient infrared regulator of the theory. 

The action ^ involves the matrices which are for- 
mally defined in the infinite Matsubara frequency space. 
In order to operate with them we have to introduce a 
cut-off for the Matsubara frequencies. Then the set of 
rules which is called J^-algebra can be established [10]. 
At the end of all calculations we tend the cut-off to in- 
finity. 

The theory ([T]) should be supplemented by the im- 
portant constraint that the combination z + T2 — 2n„r 
remains constant in the course of the RG flow. Physi- 
cally, it corresponds to the conservation of the number 
of particles [5] . In the special case of the Coulomb inter- 
action which is of the main interest for us in the paper 
the relation z + T2 — 2nvT = holds and the action ^ 
with = is invariant under a global rotation of the 
matrix Q (J^-invariance) [10]. 

The most significant physical quantities in the theory 
comprising the complete information on its low-energy 
dynamics are the physical observables a' , z'2 and z' as- 
sociated with the mean-field parameters a, Z2 and z 
of the action ([T]). The quantity a' is the conductance 
of 2DES as one can obtain from a linear response to 
electromagnetic field, n^g'j^^'^zf^/Tr is the spin suscepti- 
bility of 2DES, and z' is related with the specific heat of 
2DES [10]. Extremely important to remind that the ob- 
servable parameters a', and z' are precisely the same 
as those defined by the background field procedure [11]. 

The conductance a' is expressed in terms of the 
current-current correlations as [10] 



SUyU 



(tr[/;^,Q][/^„,g]) 



(8) 



+ 



dV((tr/;^g(r)VQ(r)trr„Q(r')VQ(r'))) 



where the limit n ^ is assumed and d denotes the 
dimension. Here and from now onwards the expecta- 
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tions are defined with respect to the theory ([T|). The 
observable is given by [5] 



rivighliBYNr dB 



(9) 



where denotes the thermodynamic potential of the 
unit volume. A natural definition of z' is obtained 
through the derivative of with respect to T [10], 



1 d n 
27rtrryA9TT' 



(10) 



To define a theory for pcrturbativc expansions we use 
the "square-root" parameterization 

' w\ 



Q = W ^hyJl-W^, W=\\ (11) 



The action ([T]) can be written as an infinite series 
in the independent fields and wl%'^^''^^^^ . 

We use the convention that Matsubara indices with 
odd subscripts ni, 77.3, . . . run over non-negative integers 
whereas those with even subscripts 712,714, .. . run over 
negative integers. 

The propagators can be written in the following form 

Rn 

_ 16^r^,,rz72 ^„,„,^C:C3^C.aj^(^^^^,^B)^.(^^^^^^B) 



2z 
a 



(27r)° 



aD'{0) - 2niReDp{0,iLUB) 



-2nlDp{0) + i2nl-l){l + j2)Dl,{0) 
+2nl{l + j2)RcD'p{0,iLUB) 



(14) 



In what follows we shall employ the dimensional reg- 
ularization scheme with d = 2 -|- 2e. Evaluating the 
momentum and frequency integrals in Eqs. (|14|) for the 
case of the Coulomb interaction {a ~ 0) we obtain for 
the physical observables p' = 4r(l - e){ATr)~'^/^ /a' , z'^ 



and z' 
1 _ 1 
P' P 



L2e r 



1 + 3^(72) - 2^254 (72) (1 + fei^Blh?)) 



Z2 



1 



(15) 



where p 
and 



ph^ 
2e 

4r(i 



l-72(2?i?,-l) + 2n272/,(cc>s//i') 
- e)(47r)-''/Vcr, fy{x) = Rc(l + ixf 



54(72) = 



1+72 
72 



ln(l +72) - 1, 



(16) 



The standard minimal subtraction scheme is the lob- 
indcpcndcnt and, therefore, is unable to treat the limits 
of vanishing [lob = Q) and strong [ujb 00) magnetic 
fields simultaneously. In order to avoid this problem 
and obtain the renormalization that correctly describes 
the two limiting cases of ojb = and ^ 00 we shall 
use the so-called cross-over renormalization scheme [12]. 

The parameter h naturally sets the momentum scale 
at which the bare parameters p, Z2 and z of the ac- 



where we introduce the notations za = 2 + r2 - 2n„r, tion IJ) are defined. The physical observables p'. 



72 = 1 + r2/z, W12 = l6nyTrTz{ni - n2)/cr, 
WB(signCi — signC2)/2 with ujb = '&nvZ2gLPBB\\/ a and 

Dp{uj,x) = [p^ + h'^ + LJ + a;]-\ 
D;{u;) = [p'' + + auj]-\ 
Dl{uj,x) = [p^ + + {I + ^2)io + x]-\ 
DpiLo) = Dpiu, 0), Dliuj) = Dliuo, 0). (13) 

The standard one-loop analysis for the physical ob- 
servables cr', Z2 s-nd z' performed at T = yields 



and z' correspond to the momentum scale h' [13, 14] 
which is determined as a'h' tv 1 = ahtr A(Q). By using 
the relation h' = h[l -\- 0{p)] we can substitute h' for h 
in Eqs. p5|) . Since the bare parameters p, Z2, and z are 
independent of h', we obtain from Eqs. p3|) the follow- 
ing one-loop cross-over RG equations in two dimensions 
(e = 0) 



a = (7 + - 
a 



d pp 
{2'kY 



dcu 



{2nl - l)j2Dliu;)Dliu;) 



+2nl-f2RcDl{Lo,tu;B)DliLj,iLOB) - (1 - a)Dl{Lj)D;{Lo) 



dp' 
drj 

drA 
drj 
dlnz' 
drj 
dfc 
drj 



= [ao(7^) + /cai(7^)]p'2, (17) 

= -p'f&o(72) + /A(72)], (18) 

co(72) + /cCi(72)], (19) 

4/c(l-/c), /c(0) = SV(l + ^')- (20) 
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Table. The function -^(72) for the cases of vanishing and strong magnetic field. Symbol lin(x) 
the polylogarithmic function. 



Er=i ^Vfc" denotes 



1+72 



+ [Anl - 1) + 72) + 2(4n2 - l)li2(-72) 



-1 

2n?, 



In 



1 - 



h 2(2n2 - 1) 
2(27i2 - 1 



li2(-72) 

1 - 



li2 



l-(2ng-l)72 



21n(l + 72), 



(2n5-l)72 



-li 



/ (2n^-l)72-l \ 
2 \^ 2nj72 J 



21ni±21. 

72 ' 



72 < 75 

72 ^ 72 



Here we have introduced 72 = Z2/Z' — 1 and variable 
rj = Inlh' where I oc corresponds to a mean-free 
path length at which, physically, the bare parameters of 
the action ([T]) are defined and 

ao(7D = -2 [1 - i^nl - 1)3,(7^)] , 60(72) = (1 + 7^)', 
CO (72) = i^nl - 1)72 - 1, «i(72) = -4^25^(72), 
61(7^) = -2n27^(l + 7^), ci(7^) = -2,127^. (21) 

The cross-over parameter fc is expressed via the quan- 
tity B = lobP. Eqs. ((iTll-llll) smoothly interpolate 
between the known results for the cases fc = and 
/c = 1 [5,6]. They were derived under assumption that 
/o' <C 1. From here onwards we omit the 'prime' sign for 
convenience. 

In general, the function /c(?7), i-c, the right hand side 
of Eq. (PO]) . is not universal. It depends on a method 
employed to derive the cross-over RG equations. The 
only universal properties of Eq. (|20p are the existence 
of two fixed points fc = and /,; = !. Fortunately, the 
qualitative features of the RG flow for Eqs. (fT7|) - ([2T|) are 
independent of the choice of fdv) provided it smoothly 
interpolates between fc = at Be^"^^ <C 1 and fc = 1 
at i3e-2i > 1. 

Let us start the analysis of Eqs. (|17p - (fT9|) from the 
case fc — 0. Then we find p{f]) oc expF(72(?7)) where 
the function ^(72) is given in Table. For all values of 
nv 72 increases monotonically with decreasing 77 and di- 
verges at 77c = T(oo; p, 72) where 

T(72;p-72) = -'e^(^^) r#i^"''^"^- (22) 

with b{u) = bo{u), p ^ p{0) and 72 = 72(0). We present 
the RG flow diagram for = 1 in Fig. H The RG 
flow is qualitatively the same for all values of Uy. As 
72 monotonous increases an initial growth of p changes 
into a decline, i.e., 2DES remains in the metallic phase 
(p 0) at large lengthscales (77 —00). 

In the opposite case /c = 1 the integration of 
Eqs. (HZll-dlSl) yields p{7j) oc exp ^"(72(77)) with the func- 
tion F(72) presented in Table. If 72 < 72 = (27i^ - 1)"^ 



(72 > 72) 72 increases (decreases) monotonous as 77 
diminishes. It reaches finally the value 72 at 77c = 
T(72;p, 72) where the function T(72;p, 72) is given by 
Eq. ([22]) with b{u) = bo{u) + bi{u). We mention that 
simultaneously p{ri) diverges at rjc- We plot the RG 
flow diagram with 7^^, = 1 in Fig. [2l Unlike the case 
fc = 0, 2DES is in the insulating phase {p 00) at 
large lengthscales for /c = 1, 72 being quenched to 72. 

In the intermediate case < /c < 1 the cross-over 
RG Eqs. p7l) - ([2T|) can be solved only numerically. If the 
parallel magnetic field is smaller than some Bx which 
is a function of p and 72 the resistance has the reen- 
trant behavior with the maximum and the minimum at 
some values of 77 (see Fig. [T|). At B > Bx the insulator- 
type behavior of the ^(77) that one expects for ;S ^ 00 
is restored. For 72 < 72 the function 72(77) has the 
maximum whereas Z2{r]) = z{r/){l + 72(7;) (spin sus- 
ceptibility) has the minimum. Their behavior does not 
change qualitatively when the parallel magnetic field 
passes through the value Bx- For 72 > 72 the max- 
imum in the function 72(77) disappears for B > Bx 
whereas 22(^7) increases monotonous for all values of B. 
As expected, the cross-over field Bx can be estimated 
as Bx ~ exp(277max) where 77niax = T(7^^''; ^,72) deter- 
mines the position of the maximum on the curve p{r]) 
at /c = 0. Here 7™'*'' is given as the root of equation 
00(72) = 0. 

p 

1 ! i 

I I 



i 1 




Fig. [2I RG fiow diagram in p versus 72 for /c = (solid 
line) and for fc = l (dashed lines) with ti,, = 1. The ar- 
rows indicate the direction towards the infrared L — > 00 
{ri —00). See text 
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We have not considered above the contribution to 
the one-loop RG equations from the particle-particle 
(Cooper) channel. It can be shown [15] that the Zeeman 
splitting due to applied parallel magnetic field does not 
affect it in the one-loop approximation. Therefore, the 
Cooper-channel contribution to the RG equations can 
be taken into account by the substitution of 00(72) — 
for 00(72) [6]. As one can check, the behavior of p{ri), 
72(7?) and 22(7/) hi this case remains qualitatively the 
same as discussed above. 

At zero temperature, h' plays the role of the inverse 
lengthscale L which is physically nothing else than a 
sample size of 2DES. If L ^ ^ D jT the temperature 
behavior of the physical observables can be found from 
the cross-over RG Eqs. (fT7 |) - (fT9| stopped at the inelastic 
length iiii. Formally, it means that one should substi- 
tute ?7T = 5 In Tl^ jD for 77 with the help of the following 
equation = 1 — [16]. The temperature depen- 

dence of the resistance in the presence of obtained 
from p(r\) in this way agrees qualitatively with the ex- 
perimental results of Ref. [7] . The detailed comparison 
between the theory and the experimental data on the re- 
sistance in the parallel magnetic field will be presented 
in the following Letter [17]. 

Also, we remind that the physical observables that we 
use throughout the paper are ensemble averaged. Even 
in a macroscopic sample of size L ^ Lin they are differ- 
ent from the measured resistance p{T) and spin suscepti- 
bility Z2 (T) due to statistically independent fluctuations 
of local conductance and electron-electron amplitudes in 
blocks of the size Lin [18]. 

Finally, we mention that the symmetry-breaking 
strain applied to 2DES with two-valleys [19] should af- 
fect the temperature dependence of the resistance in 
similar way as B|| . 

In summary, we have presented the analysis of the 
cross-over behavior of disordered interacting 2DES in 
the parallel magnetic field. The one-loop cross-over RG 
equations that smoothly interpolate between the two 
well-known limiting cases of vanishing and strong par- 
allel magnetic field allow us to explain qualitatively the 
reentrant (nonmonotonic) behavior of the resistance of 
2DES as a function of T in the presence of a relatively 
weak parallel magnetic field as well as its insulating be- 
havior at stronger fields. 
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